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We present an ab initio study of double perovskite La2VCu06, which was one of the earliest 
compounds suggested as a potential compensated half-metal. Two charge and spin configurations 
close in energy have been identified, (i) The originally envisioned spin-compensated V 4+ :d 1 /Cu 2+ :d 1 
configuration is comprised of antialigned S=l/2 cations. This state is a spin-compensated half-metal 
for moderate values of U (the on-site Coulomb repulsion strength on the metal ions) and is insulating 
for larger values of U. (ii) An unanticipated non-magnetic solution V 5+ :d°/Cu + :d 10 consists of an 
empty- and a full-shell ion, both spherically symmetric, that leads to a band insulator. This ionic 
band insulator is calculated to be the ground state at small and moderate values of U. The different 
distortions of the perovskite structure that occur for each state are central in determining the energy 
differences. Treating the Cu 2+ Jahn- Teller distortion self-consistently is particularly important for 
the magnetic solution. 



I. BACKGROUND 

Half-metals comprise an important class of ma- 
terials that has been studied extensively in recent 
years'— They are particularly interesting because 
of their anticipated applications in the field of spin- 
tronics. This type of material is conducting in one 
spin channel, but is gapped in the other. The idea 
of having such materials has been studied for nearly 
three decades£i£ This half-metallic property could 
eventually prove feasible for providing injection of 
100% spin-polarized carriers, which would provide a 
great enhancement in capabilities over current spin- 
tronics materials. 

In recent years, a twist of the argument has ap- 
peared by suggesting that one could have a half- 
metal with compensating moments, hence with the 
absence of macroscopic fields. This state is pro- 
duced with an equal number of electrons in both 
spin channels leading to no net moment, mak- 
ing the system robust against external magnetic 
fields. This situation can arise in magnetic mate- 
rials with antiferromagnetic (AF) interactions be- 
tween the magnetic ions, but such that they lead 
to partially filled bands in one spin channel and an 
insulating gap in the other. This possibility has 
been explored ab initio extensively, and a few mate- 
rials have been suggested to display this effect^— 
This property has been explored in perovskites ) 15-19 
multilayers ) 20 i 21 pnictides^ - — Heusler alloys 2 ^— 
and even diluted magnetic semiconductors. 29 Inter- 
estingly, it has been suggested that this so-called 
spin-compensated half-metals eventually could lead 
to a new type of superconductivity. 30 

One of the studied compounds has been 
La2VCuC i 6r r — where it has been suggested as a pos- 
sible spin-compensated half-metal^ but at large val- 
ues of U it has been proposed to be an insulator^ In 



this paper we explore new possible electronic phases 
in this material, apart from the previously studied 
spin-compensated solution, analyzing ab initio the 
corresponding lattice distortions that can occur and 
how they relate to each of the electronic phases con- 
sidered. 



II. COMPUTATIONAL METHODS AND 
STRUCTURE 



Our electronic structure calculations were per- 
formed within density functional theory 3 -! using the 
all-electron, full potential code wien2k22, based 
on the augmented plane wave plus local orbital 
(APW+lo) basis set. 33 The generalized gradient 
approximation^! (GGA) was used for the calcula- 
tions at U=0. To deal with strong correlation ef- 
fects we apply the LDA+U schem o 35 i 36 including an 
on-site repulsion U and Hund's coupling J for the 
3d states. We have explored various values of U for 
Cu and V, and we discuss below the different solu- 
tions that can be obtained, but our primary results 
presented here are obtained using U= 6 eV for Cu 
and U= 3.5 eV for V. For all the solutions presented, 
we have relaxed the internal positions and also the 
volume. 

La2VCuOg is studied in an ordered double per- 
ovskite structure. Due to the cation size mismatch 
between V and Cu and also due to the expected dis- 
tortions of the octahedral cage, we have optimized 
the structure (volume and internal coordinates) in 
the tetragonal space group (I4/mmm, no. 139), al- 
lowing oxygen relaxations around the cations. Two 
inequivalent oxygens occur in the structure at posi- 
tions 8h (x,x,0) and 4e (0,0,2;). 



TABLE I: Summary of interatomic distances in the two 
structures that were considered: sc (spin-compensated) 
and nm (non-magnetic) 





V-O (A) 


Cu-O (A) 


short sc 
short nm 
long sc 
long nm 


2 x 1.81 
4 x 1.83 
4 x 1.92 
2 x 1.86 


4 x 1.94 
2 x 2.01 
2 x 2.05 
4 x 2.03 



III. DISCUSSION 

We have obtained two self-consistent solutions, 
one that we will call a spin-compensated state 
(which has been discussed in the past), but also 
an unexpected non-magnetic solution. The first one 
corresponds to a formal oxidation state of Cu 2+ :d 9 , 
V 4+ :d 1 , with both atoms in a spin-half state. The 
magnetic coupling between the cations is AF, and a 
spin-compensated configuration results. The second 
is non-magnetic, with no unfilled shells, correspond- 
ing to a simple ionic Cu+:d 10 and V 5+ :d° compound. 
This emergence of high and low charge states with 
very large formal charge difference is unexpected, 
and becomes interesting as a non-magnetic competi- 
tor for the ground state. 

The differing charge and magnetic states lead to 
different types of atomic relaxation. The cation- 
oxygen distances obtained from our structure op- 
timizations are summarized in Table HI The spin- 
compensated state accommodates the V 4+ d 1 elec- 
tron in a d xy orbital by stretching the oxygen octa- 
hedron in the x — y plane, whereas the octahedron 
surrounding Cu is elongated along the z-axis to pro- 
duce one hole in the d x 2_ y 2 orbital, which becomes 
raised in energy by such distortion. Apart from re- 
laxing the internal oxygen coordinates, we have also 
optimized the unit cell volume for both states, lead- 
ing to a= 3.87A (c/a= 1.00). The internal coordi- 
nates were relaxed within LDA+U, so the particu- 
lar orbital configuration is obtained self-consistently. 
Within GGA alone, only a non-magnetic solution 
can be obtained; volume optimization in that case 
leads to a— 3.96A. 

The difference in atomic radii between Cu and V 
leads to a larger Cu-0 than V-0 distance. The ionic 
radii for Cu 2+ and Cu+ are, respectively, 0.87 A and 
0.91 A. In the case of V, the ones we need to consider 
are 0.72 A and 0.68 A for V 4+ and V 5+ , respectively. 
The elongation (respectively, shortening) along the 
c-axis of the oxygen octahedra around Cu (V) is of 
approximately 5 (6)%. For the non-magnetic solu- 
tion, the spherical symmetry of the ions leads to a 
much smaller distortion of the octahedral cages. 



The structure proposed here differs with respect 
to a recent report; 9 - where a half-metallic AF solu- 
tion with spin-half cations has been obtained, similar 
to our spin-compensated solution but with smaller 
magnetic moments. In that report a roughly undis- 
torted octahedral cage surrounding both cations was 
obtained from a structural optimization, the metal- 
oxygen distance resulting the same along the three 
axes. In particular, no Jahn- Teller distortion was 
predicted. Our results show that a structure opti- 
mization with such electronic structure distorts con- 
siderably the octahedral environment due to a Jahn- 
Teller effect, once a hole occupies one of the e g or- 
bitals of the Cu 2+ cations. According to our calcula- 
tions, the solution proposed in Ref. |9j would be un- 
stable and eventually lead, when relaxed with U= 0, 
to a non-magnetic solution. For obtaining the spin- 
compensated solution we describe here, one needs 
to use during the optimization, but it is also 
necessary to start from a Jahn- Teller distorted envi- 
ronment for both the V and Cu surrounding oxygen 
octahedra. 

Figure [T] shows the density of states (DOS) 
of the two charge configurations found. The 
spin-compensated solution (upper panel) is a half- 
metallic semimetal, with non-zero DOS at the Fermi 
level for the majority spin channel, whereas a wide 
gap appears in the minority-spin channel. For 
larger values of U, a gap opens also in the major- 
ity spin channel producing a compensated ferrimag- 
netic insulator^. However, in the reasonable range 
of U values we present here (6 eV for Cu and 3.5 
eV for V), the system is half-metallic and also spin- 
compensated. The value of the moments, both for 
Cu and V, is slightly over 0.7 /j,b inside the muffin- 
tin spheres considered, a value that is substantially 
larger than those presented in Ref. H (0.54 \ib for V 
and 0.4 \xb for Cu inside the muffin-tin spheres) and 
large enough for the S=l/2 d 1 -d 9 description to be 
realistic. 

The spin-compensated solution shows (see Fig. [TJ 
upper panel) that around the Fermi level, the V t2 g 
triplet is split by the symmetry lowering, and the 
lower d xy level becomes fully occupied. Just above 
the Fermi level lies the {d xz ,d yz } doublet. Also, 
the Cu d bands can be easily described: around the 
Fermi level there is one almost unoccupied e g band 
(of x 2 -y 2 symmetry). Its minority-spin counterpart 
is located 8 eV below the Fermi level, moved to lower 
energy by the splitting caused by U (6 eV for Cu). 

Figure [2] displays the band structures for the two 
states we are discussing. For analyzing the electronic 
structure of the spin-compensated solution, it is in- 
teresting to study only the spin-up band structure, 
since the spin-down bands shows an uninteresting 
gap around the Fermi level, as we saw above when 
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FIG. 1: Density of states (DOS) plots of the spin- 
compensated (top) and non-magnetic (bottom) solu- 
tions. Upper (lower) panels show the majority (minor- 
ity) spin density of states, in each case. A half-metallic 
solution is found for certain values of U for the spin- 
compensated solution, as shown in the top figure for U= 
6 eV for Cu and 3.5 eV for V. An insulating solution 
would occur for larger values of U (close to the unrealis- 
tic limit). The bottom figure shows the insulating, full- 
shell solution found to be the ground state non-magnetic 
solution, in this case showing the DOS for U= 6 eV for 
Cu and 3.5 eV for V. 



discussing its DOS plot. The left panel of Fig. [2] 
shows two bands crossing the Fermi level. There 
is an almost completely full V d xy band that forms 
a hole pocket around M, and also the Cu d x 2_ y 2 
crosses the Fermi level near T, being nearly unoccu- 
pied. It is clear from the band structure plot that 
a slight increase in the value of U utilized for Cu 
and/or V would make the system be an insulator, as 
was found previously 8 and described in detail above. 
The right panel of Fig. [2]shows the band structure of 



the non-magnetic solution. Just below the gap two 
occupied Cu e g bands can be identified, crystal-field 
split with respect to the Cu ta g bands, the triplet 
around -2 eV. Just above the gap, the three V t 2g 
bands can also be identified. They are crystal-field 
split with the two e g bands only slightly higher in 
energy. The gap will tend to increase for larger val- 
ues of U, but already an insulating phase occurs at 
small values of U. 

This spin-compensated solution is however not the 
ground state according to the calculated total en- 
ergy, both in the GGA (U= level) limit and also 
when the LDA+U method is utilized with moderate 
values of U (see below). The ground state is an un- 
expected non-magnetic ionic insulator, characterized 
by a d°-d 10 spherically symmetric empty shell/full 
shell combination. This non-magnetic state is in- 
sulating for a broad range of values of U, includ- 
ing the U=0 GGA limit. For the values of U we 
think are reasonable (U= 6 eV for Cu and 3.5 eV 
for V), the non-magnetic solution is more stable 
than the spin-compensated one by about 220 meV 
per metal atom. The stability of that solution gets 
greater as strong correlation effects (quantified by U 
in our calculations) are reduced. In GGA (U= 0), 
the spin-compensated solution cannot be obtained, 
as we have mentioned above. However, at larger U 
values (starting at U= 7 eV for Cu and 4 eV for 
V), the magnetic but spin-compensated solution be- 
comes insulating. It also becomes more stable than 
the non-magnetic solution, which is always insulat- 
ing. Thus, for the whole U spectrum, an insulating 
solution is always the most favorable energetically, 
ruling out the possibility of this system being a can- 
didate for showing half-met allicity in the absence of 
a total magnetic moment. 

We have also tried a different symmetry that could 
lead to new solutions^ a layered structure with al- 
ternating V and Cu layers perpendicular to the z- 
axis. This structure leads to ferromagnetic coupling 
within the V plane, produced by the distortion of the 
octahedral cage around V that leads to one electron 
occupying the {d xz ,d yz } doublet. The coupling in 
the Cu plane is AF. The total energy of such solution 
is however less stable than the non-magnetic ground 
state by about 10 meV/metal for U= 0, increasing 
for larger values of U. 

One may wonder about the contributions of the 
differing formal valences to the total energy balance, 
for the two states that compete for the ground state 
in this compound. The more highly charge differ- 
entiated V 5+ -Cu + state is lower in energy than the 
V 4+ -Cu 2+ state (except for large U), whereas a clas- 
sical picture would lead to a considerably higher 
electrostatic energy that could of course be com- 
pensated somewhat by other contributions. We just 



3 




FIG. 2: Band structures of the two states that are analyzed in the text. On the left side, the majority spin of the 
spin-compensated solution is shown. Two bands, a nearly filled V d xy valence band and an almost unoccupied Cu 
d x 2_ y 2 conduction band, overlap, leading to a half- metallic semimetal with very small Fermi surfaces. On the right 
side, the majority spin (both channels are equivalent) of the non-magnetic state is pictured, consisting of filled Cu 
3d bands and empty V 3d bands. For both Cu and V the e g -t2 g splitting is evident. The flat bands at 2 eV are La 
4/ bands; O 2p states lie below -2 eV. 



want to point out, as has been noted several times 
previously^ - — that while the formal charge state of 
ions in oxides is an extremely important and useful 
concept and label, the actual charge difference can 
be small. For these two V 5+ -Cu+ and V 4+ -Cu 2+ 
configurations, the respective actual integrated 3d 
charges within the LAPW inscribed spheres are, for 
V, 2.00 vs. 2.10, and for Cu, 8.98 vs. 8.71. Most 
interesting here is the very small difference for V, 
only 0.10 electron, and the fact that judged on ac- 
tual charge, the ionic character seems to be more 
representative of a 3+ ion in both V 5+ and V 4+ for- 
mal oxidation states. Another point of mentioning 
this is that the small actual charge difference makes 
more plausible the small energy difference that we 
find. 

Another unusual aspect is that this charge has, 
from Fig. 1, substantial contribution to hybridiza- 
tion contributions at -5 to -6 eV (with respect to 
the Fermi level). The formal valence nevertheless 
describes the overall electronic structure very well 
("fatbands" plots of wavefunction character show 
that all five V bands are unoccupied for the V 5+ 
state), and the spin moments can only be under- 
stood and discussed using the formal charge states. 



The best understanding of these effects are based on 
hybridized (Wannier) orbitals^i which include the 
effects of chemical bonding. 



IV. SUMMARY 

In summary, we have identified the possible elec- 
tronic phases of La2VCu06, one of them show- 
ing spin-compensated half-metallicity. However, the 
ground state at low and moderate values of U ac- 
cording to our total energy calculations is an un- 
expected non-magnetic insulating state, that comes 
about with the unexpected combination of valen- 
cies Cu + and V 5+ . At large values of U, a spin- 
compensated solution is more stable, but showing 
semimctallic to insulating character depending on 
the strength of on-site Coulomb repulsion. We sug- 
gest that this compound would be an interesting 
double perovskite system to study experimentally, 
to elucidate which of the two phases occurs and 
whether a metallic phase exists in a narrow window 
of parameter space. 
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